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LS I +61°303 is one of the few GeV- and TeV-emitting X-ray binaries with a prominent, well-studied modulated 
radio and gamma-ray emission. Changes in its radio morphology suggested in the past the hypothesis of a 
precessing microquasar. In 2006, a set of VLBA observations performed all around the orbit and confirming 
the fast variation in morphology were not used to study the precession because the souce was explained in 
the context of the pulsar model, the alternative model for this system. However, a recent radio spectral index 
analysis over 6.7 years from Green Bank Interferometer data at 2.2 GHz and 8.3 GHz has well confirmed the 
predictions of the microquasar scenario in LS I +61°303. At the light of these results we reanalysed the set 
of VLBA observations that constitutes a unique tool to determine the precession period and render a better 
understanding of the physical mechanism behind the precession. We improved the dynamic range of the images 
by a factor of four using self-calibration, and the self-calibrated maps reveal, in six out of ten images, a double- 
sided structure. The double-sided structure has variable position angle and switches at some epochs to a 
one-sided structure. These variations indicate a scenario where the precessing jet, inducing variable Doppler 
boosting, points close to our line of sight - a microblazar, the galactic version of the extra-galactic blazars. High 
energy observations of LS I +61°303 are consistent with the microblazar nature of this object. Moreover, we 
suggest in LS I +61°303 the first case of core shift effect observed in a microquasar. Because of this effect, well 
known in AGN, the cm-core of the jet is rather displaced from the system center. In LS I +61°303, the cm-core 
of the jet traces a large ellipse, 7 times larger than the orbit, in a period of about 28 d. Our hypothesis is that 
this ellipse is the cross-section of the precession cone of the jet at the distance of the 3.6 cm-core, and its period 
is the precession period. 



1. Introduction 



LS I +61°303 is a high-mass X-ray binary (HMXB) 
where the compact object travels through the dense 
equatorial wind of a rapidly rotating BO Ve star. The 
nature of the compact object is still unknown due to 
its large mass range, 1.4M© < M < 3.8 Mq, which 
implies cither a neutron star or a black hole. 

Two models have been proposed for this special 
HMXB, a strong and variable source at all wave- 
lenghts of the electromagnetic spectrum, from radio 
to TeV. One model assumes that the compact object 
is a non-accreting young pulsar whose relativistic wind 
strongly interacts with the wind of the Be star. The 
second model instead proposes a microquasar, that 
is, an accreting object whose steady jet, perpendicu- 
lar to the accretion disk, is occasionally travelled by 
shocks associated to transients. The peculiarity of 
LS I +61°303 is that, due to the eccentricity of the 
system, two episodes of a large mass accretion rate 
may occur along the orbit and consequently two tran- 
sients may occur per orbital period (26.496 d) (see 
references in [1]). 

In 1993, VLBI observations of LS I +61°303 showed 
that the radio emission had a structure of milliarc- 
second (mas) size corresponding to a few AU at the 
distance of 2.0 kpc [2j. However, the complex mor- 
phology in successive VLBI observations [1,13, H, H, 01 
made an interpretation in terms of a collimated ejec- 



tion with a constant position angle difficult. The 
radio morphology not only changes position angle, 
but sometimes it is even one-sided while at other 
times two-sided. This fact suggested the hypothesis 
of LS I +61°303 being a precessing microblazar (4j. 

A microblazar should be the galactic version of the 
extragalactic blazars, i.e., AGN with radio jets form- 
ing a small angle, 8, with respect to the observer's 
line of sight. Doppler boosting enhances the radiation 
from material that is moving towards the observer, 
and attenuates it, if it moves in the opposite direc- 
tion. For remarkable flux attenuation of the receding 
jet (i.e., attenuated to a level under the sensitivity of 
radio images) the structure will appear as a one-sided 
jet. A precession of the jet implies a variation of the 
angle and therefore variable Doppler boosting. The 
result is both, a continuous variation of the position 
angle of the radio emitting structure and of the flux 
density ratio between approaching and receding jet 

a. 

In known precessing X-ray binaries, the timescale 
for tidally forced precession of the accretion disk 
around the compact object, induced by the compan- 
ion star, lies within the range 8 — 22 times the orbital 
period [1, |T(|. In this context, the peculiarity of the 
variations of LS I +61°303 is their short timescale with 
respect to the predicted (8— 22) x 26.496 d. In fact, 
MERLIN images revealed a surprising variation of 60° 
in position angle in only one day Q. Even if quanti- 



eConf C111110 



2 



2011 Fermi & Jansky: Our Evolving Understanding of AGN, St Michaels, MD, Nov. 10-12 



tatively the relationship between position angle in the 
image and the viewing angle is not straightforward, 
a fast variation in position angle implies nevertheless 
clearly a fast variation in 8. The fast position angle 
variation measured with MERLIN was confirmed by 
VLBA observations [ll[ measuring in VLBA images 
a rotation of the inner structure of roughly 5°— 7° in 
2.5 hrs, that is again almost 60°/day. In that work 
some one-sided structures were taken as evidence for 
the cometary tail of the pulsar model. If the compact 
object is a pulsar, the interaction between its rela- 
tivistic wind and the equatorial wind of the Be star 
is predicted to create a bow-shock around the pulsar 
with a sort of cometary tail, i.e., a one-sided structure, 
extending away from the Be star [l2[ • 



However, recent analysis of the radio spectral in- 
dex over 6.7 years from Green Bank Interferometer 
data at 2.2 GHz and 8.3 GHz discovered new char- 
acteristics in the radio emission of LS I +61°303 [l|. 
This analysis has proven that in LS I +61°303 oc- 
curs the typical characteristic of microquasars of an 
optically thin outburst after an interval of optically 
thick emission. In microquasars, the so called tran- 
sient jet, associated to the large optically thin out- 
burst, is related to shocks travelling in a pre-existing 
steady jet, that is a slow moving continuous conical 
outflow with a composite flat/inverted radio spectrum 
(i.e., optically thick emission) [l3], [35[. The remark- 
able fact in LS I +61°303 is that, during the maximum 
of the long-term radio periodicity (4.6 yr) present in 
this source, the alternance, optically thick/thin emis- 
sion, occurs twice along the orbit, first around perias- 
tron and then again, shifted almost 0.3—0.4 in orbital 
phase towards apastron (i.e, (0.3— 0.4) x 26.496 d = 
(8—11) d after periastron) [l|. This agrees qualita- 
tively and quantitatively with the well known "two 
peak accretion/ejection model" , applied by several au- 
thors in LS I +61°303, predicting for large mass accre- 
tion rate, M , two events: one event around periastron 
and the second one shifted about 0.3 in orbital phase 
towards apastron [iH HI Il7l . 



2. Results 

2.1. Images: the fast changing structure 
from a double-sided morphology to a 
one-sided morphology 

The 8.4 GHz phase-referenced images (with respect 
to J0244+6228) shown in Fig. QjTop, well reproduce 
the structures in [ill ]. At their center, we show 
the orbit with semimajor axis of (0.22 — 0.23) mas 
(~ 0.7 x 10 13 cm at a distance of 2kpc), derived from 
Kepler's 3rd law for P = 26.496 ± 0.008 d with a mass 
for the B0 star of 17 M Q and a mass of the compact ob- 
ject of 1.4M© — 4M , and traced using the orbital pa- 
rameters from 20] . The orbit is not only much smaller 
than the radio emitting region, but also less than the 
interferometric beam of the VLBA at this wavelength. 
We remark this, because the zoomed orbit in Fig. 3 of 
(llj (zoomed to better show the orbital phase of each 
run) may unfortunately produce the wrong impression 
that the orbit is resolved whereas it is not. 

For each run are presented, in Fig. [TJ along with the 
phase-referenced ones also the self-calibrated images 
with natural and uniform weights respectively. The 
self-calibrated images were produced with the auto- 
matic procedure Muppet of the Caltech Difmap imag- 
ing package [2l|. Using self-calibration, the informa- 
tion on absolute flux density and absolute position 
are lost, but the removal of residual calibration er- 
rors improves the dynamic range up to a factor of 4 
for our maps. The higher dynamic range of the self- 
calibrated maps show in the image of run A a double 
component structure. The double structure in run 
A is clearly unlikely for a pulsar nebula, whereas it 
could be compatible with either the approaching and 
receding component of a transient jet, or with variable 
Doppler boosting along a one-sided twisted jet. Image 
B shows already in the phase referenced map with nat- 
ural weighting a double-sided morphology difficult to 
conciliate with the pulsar nebula. The self-calibrated 
map shows an intriguing twisted jet as expected for 
a fast precession. Runs C, F, I, and J show a one- 
sided jet or alternatively a pulsar cometary-tail but 
again runs D, E, G, and H show the double-sided jet 
structure. 



As a matter of fact, the radio spectral index data 
corroborate the microquasar model for LS I +61°303. 
Then, the open issue is what would produce the ob- 
served fast precession. The basis for such a study is 
to establish the precessional period. This important 
parameter could be derived by the re-analysis of the 
VLBA observations performed every 3 days over 30 
consecutive days, towards the minimum of the long- 
term radio periodicity [ll|. We present here prelimi- 
nary results of our analysis [l9j . 



2.2. Astrometry: all points form one 
ellipse 

The peak position in the phase-referenced maps are 
plotted in Fig. [5]-Bottom. It is evident that the peaks 
trace well an ellipse 6—7 times larger than the orbit. 
Our hypothesis is that this ellipse is the cross-section 
of the precession cone, at the distance of the 3.6 cm- 
core [y — 8.4 GHz) of the steady jet. 

In fact, changes in the electron energy distribution 
and the decay of the magnetic field along the con- 



eConf C111110 




= 0.865 





= 0.975 . 



= 0.090 



0=0.203 ; 



0=0.316 ; 



<£>oi 



- -Tv - 



"J I 








L5I+81°303 


G 21 Jul 2006 


Self— cal., tmif. 




m 


IB 








_ LSI+61°303 
Self-calitio+evl 

r 

1 ! ® H 
: | |j 


I 27 Jul 2006 . 





5 0-5 


LSI+61°303 
Self— cal., unif. 

• 


+ 

I 27 Jul 2006 

4- 


5 





I 30 Jul 3006 








Figure 1: Images of VLBA runs A-J at 8.4 GHz (3.6 cm) of LS I +61°303 The units on the axes are 
milliarcseconds (mas). For each run three maps are presented, the phase-referenced map (beam of (1.8—2.0) mas x 
(0.9—1.1) mas, shown in the bottom left corner), the self-calibrated map at the same resolution and the self-calibrated 
map with a beam of (1.4—1.6 mas x (0.6—0.9) mas. Contour levels for all maps are —4, 4, 5.66, 8, 11.3, 16, 22.6, 32, 
45.2, 64, 90.5, 128, 181 a (with la — 0.1 mjy/beam for Top and Bottom maps, la ~ 0.07 mjy/beam for Middle maps). 
At the center of the phase-referenced images is traced, in scale, the orbit (see Sect. 2). 
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ical outflow forming the steady jet, imply that the 
critical frequency, i-Wak ; of the synchrotron spectrum 
varies along the jet [22L [35| . In microquasars, the 
^reak for the part of the steady jet closest to the en- 
gine (i.e., i^brcaki in Fig- 2) lies in the infrared [23[, 
whereas in AGN the observed turnover lies in the mil- 
limeter range 22]. Observing at longer wavelengths 
than infrared, i.e., in the radio band, gives rise to 
two results: a flat spectrum and the "core shift". 
Multi-wavelength observations result in fact in the 
flat/inverted spectrum discussed in the introduction, 
typical for a steady jet. Imaging a steady jet at one ob- 
serving frequency i^obs; gives rise to the effect known as 
"core shift" [25[ , with the displacement from the cen- 
ter as a function of the observing frequency f bs- At 
z^obs the emission of the segment will dominate, whose 
spectrum peaks at that frequency plus small contribu- 
tions from neighboring segments (see Fig. 1 in (24|). 
With the engine being close to the "infrared-core" , it 
is clear that the 3.6 cm-core can be rather displaced 
from the center (like segment 2 is displaced with re- 
spect to the center in Fig. 2 Top-Left) [H [H HI]. 
If a jet is pointing towards us (i.e., a microblazar as 
in Fig. 2 Top-Right) the core will be dominated by 
the approaching jet, because of Doppler boosting. If 
the jet is precessing, then the core will describe an el- 
lipse (see both ellipses, the observed ellipse in Fig. 2- 
Bottom and the predicted ellipse for the blueshifted 
segment-2 in Fig. 2 Top- Right). A transient jet, on 
the contrary, can have any shift from the center, de- 
pending on the velocity and the time elapsed from 
the transient. The peak of map E, even if displaced 
from the other ones, and therefore likely affected by 
the approaching component of the transient jet, is at 
a position angle which is consistent with the other 
peaks. 

Note as after 27 d, the peak of run I is rather close 
to completing the cycle, i.e. to overlap with the peak 
of the starting run A. The same occurs for peak of 
run J, that 27 d after run B, nearly overlaps with 
peak B. The peak of run J, 30 d after run A, is clearly 
rather displaced from the A peak. The period seems 
therefore to be in the range (27—28) d. 



3. Conclusions and Discussion 

LS I +61°303 is one of the few established massive 
X-ray binaries that emit in the high and very high 
energy range. It is formed by a compact object of un- 
known nature (black hole or neutron star) travelling 
with a period of 26.496 d around a Be star. Two mod- 
els have been proposed in the past: a millisecond pul- 
sar and a microquasar model. Recent analysis of the 
radio spectral index has proven in LS I +61°303 the 
typical characteristic of microquasars of an optically 
thin outburst after an interval of optically thick emis- 
sion twice along the orbit, as predicted by the two 
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Figure 2: Top-Left: Superposition of individual spectra, 
each with a different break frequency associated with 
different segments of a steady jet. Top-Right: For a 
precessing microblazar, i.e., for small angle of ejection 
with respect to the line of sight, the core component, 
dominated by the approaching jet contribution because 
of Doppler boosting, will describe an ellipse during 
precession. Bottom: Astrometry of consecutive peaks of 
VLBA 8.4 GHz maps for runs A-J with the orbit drawn 
in scale at an arbitrary distance. The quality of the 
astrometry is determined from the scatter in positions of 
the check source J0239+6005 [l]J. The peak of map E, 
even if displaced from the other ones, and therefore likely 
affected by the approaching component of the transient 
jet, is at a position angle which is consistent with the 
other peaks. The fact that the precessing compact object 
is moving in an orbit with semimajor axis of ~0.2 mas 
introduces additional variations [19j . 



peak microquasar model. 

The results of our re-analysis of VLBA observa- 
tions of LS I +61°303 presented here are that the 
radio emission has in several images a double-sided 
structure (see. Fig. 1). The astrometry show that 
the peaks of the images trace a well defined ellipse 
in (27—28) d (see Fig. 2). The pulsar model explains 
neither the double-sided morphology nor the observed 



eConf C111110 



2011 Fermi & Jansky: Our Evolving Understanding of AGN, St Michaels, MD, Nov. 10-12 



5 



change from double sided to a one-sided structure. 
The microquasar model can explain them with vari- 
able Doppler boosting, i.e., with a precessing jet. The 
cm-core of a precessing steady jet pointing close to 
our line of sight, as in a microblazar, is expected to 
describe an ellipse during the precession. During the 
transient jet phase there will be an additional shift due 
to the approaching jet component. We conclude there- 
fore that the precession period is the time of (27—28) d 
necessary to complete the ellipse. 

A precession period of (27—28) d for the accretion 
disk in LS I +61°303 induced by tidal forces of the Be 
star (loj would require the unrealistic value for the size 
of the accretion disk of 0.5 — 0.8 x 10 13 cm, i.e., nearly 
the semi-major axis, and can therefore be ruled out. 
On the contrary, for a slow rotating compact object, 
a precession period of 28 d induced by Lense-Thirring 
effect (i.e., frame dragging produced by the rotation 
of the compact object) could be compatible (loj . 

High energy observations of LS I +61°303 are con- 
sistent with the microblazar nature of this object. Si- 
multaneous X-ray and VHE observations during an 
outburst of LS I +61°303 resulted in a correlation, 
indicating a simultaneity in the emission processes 
p6^ . In particular, with respect to the comparison 
with blazars, it has been noticed [l(| as the two fluxes 
result in -Fvhe cx with 77 = 0.99 in agreement with 
the correlation observed in blazars 27]. In addition, 
the respective values of the photon indices seem to 
be comparable. During a bright flaring event of the 
blazar 1ES 2344+514, VERITAS measured a photon 
index of L = 2.43± 0.22, whereas for the X-ray emis- 
sion, RXTE and Swift-XRT measured a hard pho- 
ton index of -1.9 [28]. For LS I +61°303: MAGIC 
had a F = 2.7 and the X-ray emission detected by 
XMM-Newton and Swift-XRT a harder photon index 
of 1.5-1.8 [26|. In blazars the X-ray emission is due 
to synchrotron, and VHE is synchrotron self-Compton 
(SSC) [29(. As discussed in [30(, several authors have 
in fact explained the X-ray excess around apastron in 
LS I +61°303 with synchrotron and the VHE with ei- 
ther external inverse Compton (EIC) or SSC emission 
[3TI HH. In particular, the X-ray /VHE correlation 
in LS I +61°303 [26] is compatible with a one zone 
leptonic particle population producing the emission 
[2a . l32j . Concerning high energy emission in the GeV 
range, detected e.g., with Fermz-LAT, GeV emission 
is seen all along the orbit. In fact, electrons from 
the steady jet can always upscatter stellar UV pho- 
tons to GeV energies (i.e., EIC see [Hj]). Nevertheless, 
more energetic particles from the transient jet could 
in addition also produce GeV emission via EIC and 
SSC. Intriguingly, the spectrum measured by Fermi- 
LAT shows, in addition to a power law with a cut-off 
around 6 GeV, upper limits possibly compatible with 
the spectrum measured with MAGIC and VERITAS 
(see e.g., Figs. 2 and 3 in [33| and discussion in [30j]). 
As a matter of fact, there is an interesting increase in 



the overall flux level observed with Fermi-LAT after 
March 2009 together with a broadening of the peak 
shape [33(1 . The up to now observed Fermi-LAT vari- 
ations are therefore consistent with a long term varia- 
tion. Similarly, strong variations are observed at very 
high energies. With VERITAS the source went from 
being detected around apastron to becoming quiescent 
between 2008 and 2010 [Hj]. With respect to the long- 
term, 4.6 yr, radio periodicity mentioned in Sect. 1, 
the insufficient temporal coverage at high energy, ev- 
ident in Fig. 1 of [20], does not allow at the moment 
a more closer comparison. However, a different trend 
with respect to the 4.6 yr radio periodicity is expected. 
In fact, the high energy emission is indeed related to 
the steady jet but its peak seems to occur during the 
transient jet and not during the steady one (Fig. 4 
in (35[). As a matter of fact the steady jet, related 
to the radio periodicity, is not always followed by the 
transient one. A timing analysis of the transient jet is 
in progress. 
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